Optical properties of silicon nanoparticles in the presence of water: A first principles 

theoretical analysis 
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We investigate the impact of water, a polar solvent, on the optical absorption of prototypical silicon 
clusters with oxygen passivation. We approach this complex problem by assessing the contributions 
of three factors: chemical reactivity; thermal equilibration and dielectric screening. We find that the 
silanone (Si=0) functional group is not chemically stable in the presence of water and exclude this 
as a source of significant red shift in absorption in aqueous environments. We perform first principles 
molecular dynamics simulations of the solvation of an oxygenated silicon cluster with explicit water 
molecules at 300 K. We find a systematic 0.7 eV red shift in the absorption gap of this cluster, 
which we attribute to thermal strain of the molecular structure. Surprisingly, we find no observable 
screening impact of the solvent, in contrast with consistent blue shifts observed for similarly sized 
organic molecules in polar solvents. The predicted red shift is expected to be significantly smaller for 
larger Si quantum dots produced experimentally, guaranteeing that their vacuum optical properties 
are preserved even in aqueous environments. 
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I. INTRODUCTION 

Molecular sensing technology currently relies on fluo- 
rescent organic molecules as optical tags. However, these 
organic dyes have some deficiencies. They absorb light 
only at specific frequencies, making them inefficient emit- 
ters. Their optical emission cannot be tuned readily, and 
many molecules with unrelated structures and chemistry 
are required to produce a range of emission wavelengths. 
Furthermore, prolonged irradiation can lead to chemical 
or structural decomposition, altering or even destroying 
desirable optical properties during experiments. These 
deficiencies may be overcome by using semiconductor 
quantum dots (QDs) to develop new optical sensing tech- 
nology. 

QDs arc pieces of bulk semiconductor in the 
nanometer-size regime. Their small size bestows them 
with optical properties that lie somewhere between those 
of molecules and those of the corresponding bulk ma- 
terials. Like their bulk counterparts, QDs have robust 
atomic structures. They absorb light efficiently over a 
wide range of wavelengths, while emitting light of a char- 
acteristic frequency. In addition, due to quantum con- 
finement effects, the optical properties of QDs may be 
tuned by controlling the size of each nanoparticle. These 
useful properties make QDs highly competitive with or- 
ganic dyes for use as optical tags. 

QDs made of silicon (Si) are highly desirable, given the 
abundance of this element and its use in the electronics 
industry. Also, in contrast with existing QD materials 
such as cadmium selenide (CdSe), porous silicon (one 
of the precursors of Si QDs) is both bio-stable and non- 
toxic, jj Currently, the use of CdSe QDs as optical tags in 
biological environments requires coating with many lay- 
ers of materials, to increase bio-stability, and shield the 
biological environment from their inherent toxicity. Si 
QDs do not require any coating and therefore the over- 
all size of Si QDs is typically much smaller than that 
of coated CdSe QDs. This should permit the use of Si 



QDs, which possess desirable optical properties and may 
be applied in aqueous environments, 2] to probe bio- 
chemical processes currently inaccessible to CdSe QDs 
- specifically processes which involve diffusion through 
intra-cellular membranes. 

However, at present, Si QDs of a particular size, syn- 
thesized using different approaches, have a large dis- 
tribution of measured optical properties, e.g. photolu- 
minescence from UV to red. |3( To shed light on the 
discrepancies between different experiments, first prin- 
ciples calculations have illustrated the imp ortant role 
of surface chemistry U, IE IE IE IE Il0| and struc- 



ture [il El EE H El on the optical properties of Si 
QDs. From this body of work, it is now clear that the 
optical properties of Si QDs are not determined by size 
alone. 

In this work we analyze the impact of water on the 
optical absorption properties of Si nanoparticles. In the 
synthesis process or the application environment, water 
may be present, in vapor or liquid form, as a solvent or 
a contaminant. Furthermore, in applications to sensing 
in biological environments, water would be ubiquitous. 
Given the known impacts of water - a polar solvent - on 
the absorption properties of solvated organic molecules, 
and the absence of equivalent information for inorganic 
solutes, analysis for Si nanoparticles in water is clearly 
necessary and important. 



II. CONSEQUENCES OF SOLUTE-SOLVENT 
INTERACTION FOR OPTICAL PROPERTIES 

In our study, we consider the solvation of an optically- 
active species in a polar solvent. The strength of the 
solute-solvent interaction determines how much the op- 
tical properties of the total solvated system deviate from 
those of its constituents: the isolated solute and the sol- 
vent. If the electronic states involved in optical processes 
can be clearly associated with the solute (rather than 
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with the solvent) then we may speak of the impact of 
the solvent on the solute, and that is the focus of this 
paper. j^jj We consider three solute-solvent interactions 
contributing to changes in the optical properties of the 
solute: chemical reactivity; thermal equilibration; and 
dielectric screening. 

Dielectric screening refers to the impact on the molec- 
ular and electronic structure of the solute due to screen- 
ing of Coulomb interactions caused by the finite electri- 
cal susceptibilities of the solvent and solute. From the 
electrostatic point of view, this is a self-consistent phe- 
nomenon, where the charge density of the solute estab- 
lishes a polarization of the solvent, which in turn im- 
pacts the solute charge density until some equilibrium is 
reached. Concomitantly, long-range electron correlation 
leads to attractive dispersion forces between the solute 
and solvent charge densities. The impact of these screen- 
ing contributions on optical properties is explained in 
detail in Sec. 11111 and the impact of screening on the Si 
clusters studied here is discussed in Sec. lVIIll Sections lTVI 
and contain details of our first principles calculations 
and their application to Si clusters in vacuo, respectively. 

An investigation of the reactivity of various Si clus- 
ters with water is given in Sec. IVII Clearly, if the solute 
is chemically reactive in the presence of the solvent, the 
resulting reaction product may have a different molecu- 
lar structure or composition, and, consequently, a com- 
pletely different optical signature. 

Finally, in Sec. IVIll we examine the effect of strain in- 
duced by thermal fluctuations on optical properties of Si 
clusters. Thermal equilibration with a reservoir of sol- 
vent molecules will activate thermally accessible vibra- 
tional modes of the solute. Already, for Si clusters in 
vacuo at zero temperature, there is clear evidence that 
strain in the molecular structure, e.g. due to surface re- 
laxation, has a large im pact on the optical properties of 
the system. 
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III. SOLVENT SCREENING AND OPTICAL 
ABSORPTION 

In order to study the impact of water on the opti- 
cal absorption of Si nanoparticles, it suffices to assess 
the impact of solute-solvent interactions on the electronic 
ground state and on the first excited state. The solva- 
tion energy |4pj may be different for ground and excited 
states, and this is the origin of so-called solvation shifts 
(Fig. 2]) . If the solvation energy of the excited state is 
larger (smaller) than that of the ground state, we see a 
red (blue) shift in the absorption energy when compared 
with the unsolvated system. 

One source of this difference in solvation energies is di- 
electric screening. For a large class of systems, probabil- 
ity distributions of excited state wave functions are more 
diffuse than those of the ground state, and consequently 
the system has a higher polarizability in the excited state 
than in the ground state. In general, the more polariz- 
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FIG. 1: Schematic energy level diagrams for the electronic 
states involved in optical absorption for vacuum and solvated 
systems. Absorption energies required for transitions from 
ground to excited states are indicated with arrows. The sol- 
vation energy causes shifts in the total energy of the ground 
and excited states. The difference in solvation energies of the 
ground and excited states leads to a red or blue shift in the 
absorption energy with respect to the absorption energy in 
vacuo. 



able a system, the larger its solvation energy. Therefore, 
if the excited state is screened by the solvent, this larger 
solvation energy will lead to a red shift in absorption. 
This is often the case for non-polar solvents. 

In polar solvents a significant contribution to the 
screening comes from the alignment of molecular dipoles 
and occurs on the timescale of molecular motion. Upon 
absorption of a photon, there is not sufficient time for this 
orientational polarization to screen the resulting change 
in the electronic charge density of the solute. Hence, the 
resulting non- equilibrium excited state has a reduced sol- 
vation energy with respect to the ground state, leading 
to an overall blue shift in absorption (Fig. [TJ. 

In practice, both the fast and slow components of the 
screening response of the solvent are present, leading to a 
net solvent shift. Water is a highly polar solvent, with an 
average dipole moment per molecule of ~ 1.9 Debye in 
the gas phase. The electrical susceptibility of water can 
be separated into a static component \o and an instanta- 
neous component Xoo • (xo = e o ~ 1 , where eo is the static 
dielectric constant, and Xoo = £oo — 1 = n> — 1, where 
is the optical dielectric constant and n is the refractive 
index). For liquid water at room temperature, Xoo = 0.78 
and xo ~ 77. Experimental estimates of the Debye re- 
laxation time for water provide an estimate of the time 
taken for reorganization of the orientational polarization 
in response to an external electric field. Current theories 
and experiments agree that there are two Debye relax- 
ation times for water, one fast and one slow. The fast re- 
laxation time has been tentatively associated with single 
molecular reorientation and has values from 0.2 - 1.5 ps 
at ambient conditions. |!6lll7| Therefore we can assume 
that electronic absorption processes (occuring on 0.1 fs 
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FIG. 2: Schematic energy level diagrams of the electronic 
states considered in this paper, where single particle orbitals 
are indicated by horizontal lines, energy increases in the ver- 
tical direction, and up and down spin channels are separated 
with a dashed line, (a) The spin degenerate ground state, 
indicating the highest occupied molecular orbital (HOMO) 
and the lowest unoccupied molecular orbital (LUMO). (b) 
The first excited state with triplet symmetry in a single ref- 
erence picture, where the highest occupied (lowest unoccu- 
pied) molecular orbital of this excited state configuration, la- 
beled HOMO* (LUMO*) can be associated with the LUMO 
(HOMO) of the ground state, (c) the first excited state 
used to approximate singlet symmetry, with the HOMO* and 
LUMO* in the same spin channel. 



time scales) are instantaneous and that the excited states 
correspond to a non-equilibrium nuclear configuration. 
In addition, the large difference in magnitudes between 
Xoc and xa increases the energetic cost of the excited 
state and generally leads to blue shifts in the absorption 
spectra of optically-active polar organic molecules, e.g. 
acetone (0.2 eV) and pyridazine (0.5 eV). In this 
paper, we will see (surprisingly) that these characteristic 
blue shifts are absent (or at least negligible in the pres- 
ence of thermal strain) for small inorganic Si clusters (see 
Sec. EDI. 



IV. FIRST PRINCIPLES CALCULATIONS 

The optical properties of the various Si clusters ex- 
amined in this paper are known to be particularly sen- 
sitive to their molecular structure. |13| Furthermore, 
the properties of the solvent and the nature of the 
solute-solvent interaction can be quite complex, involv- 
ing charge-transfer, hydrogen-bonding and long-range 
screening. To capture all of this in our theoretical ap- 
proach, we need highly accurate estimation of the elec- 
tronic structure of these systems. For these reasons we 
adopt first principles computational techniques. 

We make use of density functional theory (DFT) un- 
der the generalized gradient approximation (GGA). To 
simulate the impact of finite temperature and explicit 
solute-solvent molecular interactions, we use first prin- 
ciples molecular dynamics (FPMD), employing the Car- 
Parrinello method, (2^ treating the nuclei in the system 
as classical particles acted upon by Coulombic inter-ionic 
forces and quantum forces derived from the electronic 
structure. 



The majority of electronic structure calculations in 
this paper were performed using the plane-wave pseudo- 
potential codes GP [2l| and abinit. [23 We also used 
the all electron code GAUSSIAN 98. |23| FPMD sim- 
ulations were carried out using GP. The DFT/GGA 
plane- wave pseudo-potential calculations make use of the 
PBE exchange correlation functional [24[ and use a basis 
of planewaves consistent with periodic boundary condi- 
tions on a simulation box of a given size. All calcula- 
tions are performed for a single point in the Brillouin 
zone, namely, the r-point. In these plane- wave calcula- 
tions, only valence electrons are included, with atomic 
cores represented by nonlocal, norm-conserving pseudo- 
potentials. [2^, |2(J The Gaussian-basis all electron calcu- 
lations make use of the Becke functional for exchange 
and the PW91 correlation functional. |28j These calcula- 
tions are performed for finite systems using a 6-31 1G** 
Gaussian basis set. 

A plane-wave kinetic energy cut-off of 70 Ry 41] pro- 
duces converged results for the electronic structure of the 
oxygen-containing Si clusters considered here and for wa- 
ter in the condensed phase. For calculations involving Si 
clusters in vacuo, simulation cell dimensions are chosen 
appropriately to reduce the impact of finite size effects on 
the results. This is particularly important in calculations 
on small molecules where the spatial extent of unoccu- 
pied electronic states can be much larger than the atomic 
dimensions of the molecule. 

As discussed in Sec. 11111 to determine the absorption 
properties of a given system, we require explicit knowl- 
edge of the equilibrium electronic ground state and the 
non-equilibrium excited state. In practice, the absorp- 
tion spectrum of a physical system is composed of many 
allowed electronic transitions and associated vibrational 
transitions. In this paper, we concentrate only on the 
lowest possible electronic transition in the absence of vi- 
brational excitations. This marks the onset of optical 
absorption and is equivalent to the band gap of a solid 
state system. We shall frequently refer to this energy as 
the absorption gap, E g . In assessing the impact of water 
on the optical properties of Si clusters, we assume that 
a change in the absorption gap is indicative of an overall 
shift in the absorption spectrum of these systems. 

Estimation of the absorption gap using first princi- 
ples methods can be difficult. Electronic excitations are 
at heart many-body processes, and so, are difficult to 
describe using effective single-particle DFT orbitals, es- 
pecially since most correlation functionals are designed 
specifically for the electronic ground state. One com- 
mon DFT estimate of E g is the difference between the 
Kohn-Sham eigenvalues of the highest occupied molecu- 
lar orbital (HOMO) and the lowest unoccupied molecular 
orbital (LUMO) of the ground state electronic structure 
[Fig- 121 (a)]- For finite systems in vacuo, This HOMO- 
LUMO gap has been shown, to follow the same trends, 
as a function of size, as more accurate estimates of the 
absorption gap computed using quantum Monte Carlo 
(QMC) calculations. [29J 



However, in determining the impact of dielectric 
screening on the optical properties of a solvated sys- 
tem, we note that the DFT HOMO-LUMO gap is not 
an appropriate estimator of the absorption gap. Since 
the LUMO is never occupied, the HOMO-LUMO gap 
does not include physical information on the change in 
solute charge density and corresponding response in sol- 
vent polarization upon electronic excitation. For solvated 
systems, we choose to employ the ASCF approach, to 
determine the absorption gap. Here, we define E g — 
Et ot (X ) - E tot (X ), where E^ ot (X ) is the total excited 
state energy and E to t(Xo) is the total ground state en- 
ergy, both evaluated for the ground state molecular struc- 
ture Xq. By construction, E% ot includes the impact of 
electron density variation associated with the solute ex- 
cited state and any screening impacts from the solvent in 
response to the solute. The importance of the ASCF es- 
timator of Eg for solvated systems is discussed in Sec. lIXI 

For the systems considered here, the first allowed elec- 
tronic transition from the ground state is to a state with 
singlet symmetry. The representation of this singlet state 
requires more than one reference electron configuration. 
Since DFT is an effective single particle method we ap- 
proximate this excitation with the single reference state 
described in Fig. [2^c). This approximate "singlet" does 
not have the same symmetry as the true singlet, but it 
does contain a realistic Coulomb interaction between the 
states near the Fermi level. Furthermore, a QMC calcu- 
lation of the absorption gap of the silane molecule (S1H4) 
using this single reference approximation has been shown 
to be accurate to within 0.1 eV. |30j 

However, it is well known that numerical convergence 
of the DFT "singlet" first excited state energy is inher- 
ently difficult. A more practicable energy to calculate 
is that of the first excited electronic state with triplet 
symmetry [Fig [3(b)]. This triplet state has a reduced 
Coulomb energy when compared with the "singlet" , due 
to the exchange repulsion that keeps like-spin electrons 
apart. For reasons of efficiency and robustness, we choose 
to compute the ASCF estimate of E g using the triplet 
excited state. In practice, using either the "singlet" or 
triplet excited state energy does not qualitatively change 
our results, and we provide evidence for this in Sec. IIXI 
by comparing with some "singlet" state calculations. In 
the rest of this paper, we shall refer to the "singlet" state 
without quotation marks. 



V. SILICON CLUSTERS IN THE ABSENCE OF 
WATER 

The Si clusters studied in this paper are presented in 
Fig. [3 and the corresponding estimates of the absorption 
gaps in vacuo are presented in Table [I] Despite definite 
inadequacy in the DFT estimates of the absorption gap, 
when compared with accurate QMC calculations, there 
are still consistent trends in comparison of the gaps of 
Sis clusters with different passivants. This has also been 
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FIG. 3: The structures of various silicon clusters and cor- 
responding probability distributions (green) of the highest 
occupied molecular orbital (HOMO) and lowest unoccupied 
molecular orbital (LUMO). All electron density isosurfaces 
contain 30% of the total charge of the given state. Si atoms 
are colored blue, O red and H white. 



demonstrated in previous work over a wide range of Si 
cluster sizes. Q Also, the smaller systematic differences 
between the two sets of DFT gaps are due to the choice of 
exchange correlation functional and whether the pseudo- 
potential approximation is adopted. We exploit the con- 
sistent trends in these calculated gaps, which allow quali- 
tative comparison of the gaps of different systems as long 
as they are computed using the same level of theory. 

Our reference structure is the "pure" hydrogen- 
saturated Si5Hi2 [Fig. dfa)]- We also consider various 
types of surface passivation containing oxygen. The first, 
Si5HioO(I) contains a silanone (Si=0) group [Fig.|3Jb)]. 
The departure from the local tetrahedral symmetry of 
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reactants transition state 



product 



TABLE I: Absorption gaps (eV) and dipole moments (Debye) 
for various silicon clusters. Gaps estimated using: HOMO- 
LUMO DFT/PBE pseudo-potential calculations (PBE-pp); 
HOMO-LUMO DFT/PW91 Gaussian basis set, all-electron 
calculations (PW91-ae); and "singlet" ASCF QMC calcula- 
tions (QMC). Dipole moments estimated from PBE-pp charge 
densities (d vac )- Also shown are dipole moments calculated 
in the presence of a continuum dielectric solvation model for 
water (d so i v ) and the angular change in direction (8) from the 
unsolvated dipole moment. (The symmetry of SisHi2 guar- 
antees that its dipole moment is zero.) 
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each Si atom to a planar sp 2 symmetry for the SiHO 
portion of the molecule produces a localized HOMO and 
LUMO on this functional group. The localization effec- 
tively creates a defect state within the "band gap" of the 
molecule reducing the absorption gap by almost 50%. We 
refer to previous work for a discussion of this effect for 
Si clusters of various sizes. 

We next consider two other forms of oxygen passivation 
which maintain, to a greater degree, the local tetrahedral 
symmetry of the Si atoms in the cluster. Si5HioO(II) 
contains a bridging oxygen atom between neighboring Si 
atoms at the surface [Fig. [3c)]. This causes some strain 
in the angle between the surface Si atoms bonded to the O 
and the core Si atom, reducing the tetrahedral angle from 
the ideal value (109.5°) to 70.4°. However, the preserva- 
tion of approximate tetrahedral symmetry around each 
Si atom prevents the same degree of localization of the 
HOMO and LUMO as we see in Si 5 Hi O(I). Conse- 
quently, the impact on the absorption gap is diminished 
and, for Si 5 Hi O(II), E g is 25% less than that of Si 5 Hi 2 . 
SisHnOH contains a hydroxyl (OH) passivant, which has 
a reduced impact on the molecular structure in terms of 
strain, in comparison with Si5HioO(I). However, there is 
still a reduction in the absorption gap due to the electro- 
negativity of the O atom, which attempts to localize the 
states near the Fermi level around it [Fig. Eld)]. 

In a first approximation to the impact of water on the 
structure of these Sis clusters, we use a continuum solva- 
tion model [3TI [3^ to simulate the electrostatic impact 
of the solvent. We find almost no change in the atomic 
structure of these molecules in the presence of this solva- 
tion model. The most noticeable change is a 1% increase 
in the bond length of the silanone group in Si5HioO(I). 
However, an examination of the impact of this solvation 
model on the dipole moments of each molecule indicates 
a consistent tendency to increase the dipole moment (Ta- 
ble P) . Given that the atomic structure remains unper- 
turbed in the presence of the solvent polarization, this in- 
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FIG. 4: The reaction paths of Si 5 Hi O(I) (red) and 
Si5HioO(II) (black) with one water molecule, indicating the 
activation energy E ac ^ and the overall energy of reaction 



TABLE II: Variations in the HOMO-LUMO gap (eV) during 
the reaction with one water molecule for SisHioO(I) and (II). 
Calculated using PW91-ae. 
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Si 5 H 10 O (I) 
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4.0 


5.3 


SigHioO (II) 


4.7 


4.1 


4.7 



dicates that the increase in the dipole moment is caused 
by polarization of the electron charge density. 



VI. REACTIVITY OF OXYGENATED SILICON 
CLUSTERS WITH WATER 

The most striking impact of the solvent on a solute 
would be a chemical reaction producing a more stable 
chemical species, which we could expect to have quite 
different optical properties to its precursor. Of the four 
clusters considered in Section SisHioO(I) readily re- 
acts with water to produce a dihydroxide according to 
the reaction: 



Si 5 Hi O(I) + H 2 => Si 5 H 10 (OH) 2 



(1) 



Using GAUSSIAN 98, we determined a 0.1 eV activa- 
tion energy for this reaction. This involved finding the 
total energy of the transition state for the reaction, by 
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structural optimization in the direction of the eigenvec- 
tor of the Hessian matrix with lowest negative eigenvalue. 
This reaction profile is compared, in Fig. 01 with the 
same reaction for the stoichiometrically equivalent bridg- 
ing structure Si5HioO(ll), which has a larger activation 
energy of 0.9 eV. This analysis indicates that the survival 
of silanone groups in water is highly unlikely, given the 
comparable size of the activation energy with the ther- 
mal energy at 300 K (i.e., ~ 0.03 eV). Also, this reaction 
is exothermic, and the energy released (~ 1.8eV) would 
fuel the breakdown of other silanone groups by water. 

Previous work by Zhou and Head [33j | indicates a larger 
activation energy of 0.3 eV for the silanone group's reac- 
tion with water. However, the level of theory (MP2) and 
smaller system size (one Si atom) in that work may lead 
to a larger activation energy than the one found here. 
Furthermore, the 1% increase in the Si=0 bondlength 
in the presence of a polarizable solvent model (Sec. [Vj) 
indicates that in liquid water this activation energy may 
be further reduced due to polarization of the Si5HioO(I) 
molecule and weakening of the Si=0 bond. A recent 
study of the reactivity of several silica clusters with wa- 
ter by Laurence and Hillier [34j, using the B3LYP func- 
tional and a comparable basis to ours, reports activation 
energies of ~ 0.9 eV for the reaction of the Si-O-Si group 
with water, in good agreement with our result for the 
reaction of Si5HioO(II) with a water molecule. 

The reaction of Si5HioO(I) with water leads to a prod- 
uct with a larger absorption gap. The restoration of 
the local sp 3 symmetry around the Si atom of the Si=0 
group relaxes the strain on the electronic structure, and 
the localized states present in the reactant disappear. 
This increases the DFT HOMO-LUMO gap from 3.8 eV 
(in the presence of a water molecule) to 5.3 eV, compa- 
rable to that of the single hydroxide. 

We have also observed this reaction to take place in 
the course of a FPMD simulation of Si5HioO(I) in a sim- 
ulation cell containing 57 water molecules. This reaction 
occurred while the system temperature was being raised 
to 300 K, within the first 0.2 ps of the simulation. There- 
fore, we exclude this cluster from further tests on absorp- 
tion in the presence of water. Given that the product of 
this reaction is a hydroxide, which is significantly more 
stable than the silanone, we regard the hydroxide cluster 
as a chemically stable species in water. 

The absorption gap of Si5HioO(II) remains essentially 
the same at 4.7 eV, following the reaction with water. We 
regard this bridged cluster as chemically stable given its 
larger activation energy (0.9 eV), and smaller exothermic 
energy (0.6 eV) in comparison with that of the silanone 
cluster. Furthermore, as we see in Sec. IVIII during a 
4.5 ps FPMD simulation we do not observe the hydrox- 
ylation of this molecule when solvated in liquid water, 
despite large fluctuations in the cluster's average kinetic 
energy which were equivalent to temperatures exceeding 
1000 K [Fig. EJa)]. 
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FIG. 5: (a) The system temperature Ttot (black) and the 
instantaneous average kinetic energy (converted to units of 
temperature) of the Si atoms Tsi as a function of time dur- 
ing an FPMD simulation of Si 5 Hi O(II) in water, (b) The 
absorption gap of the solvated system (Si cluster and water 
molecules) computed at the indicated times during the course 
of the FPMD simulation (black) and the corresponding gap 
of the same Si cluster structure without the water molecules 
(red). Also shown is the gap of the Si cluster in its relaxed, 
zero temperature structure (green). The arrow indicates the 
particular snapshot used to produce the electronic structure 
information in Figs. and |7|and Tabic ITTT1 



VII. IMPACT OF FINITE TEMPERATURE ON 
OPTICAL PROPERTIES OF SILICON CLUSTERS 

We simulate the solvation of the molecule SisHioO(II) 
in water using FPMD. The simulation cell dimensions are 
(12.169 A x 12.169 A x 12.169 A) and the cell contains 
57 water molecules in addition to the solute. These cell 
dimensions were optimized to produce ambient pressure 
conditions using a classical MD simulation. The final 
configuration of this classical simulation was used as the 
starting point of the FPMD run. 

The initial 0.17 ps of the FPMD simulation evolve un- 
der the influence of a velocity rescaling thermostat after 
which the temperature of the system reaches 300 K. At 
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this point the thermostat is removed and the simulation 
is allowed to evolve with constant total energy. 

We notice that the molecular structure of the solute is 
greatly perturbed during this simulation. This perturba- 
tion is caused by the non-equilibrium initial configuration 
that is used to start the simulation. We assume that the 
classical forces used to generate this initial configuration 
have some mismatch with the quantum mechanical forces 
present in the FPMD simulation. This mismatch greatly 
increases the average kinetic energy of the solute, as can 
be seen in Fig. EJa) (the average kinetic energy of the 
Si atoms in the solute, T$i, is shown here converted to 
temperature units for comparison with the system tem- 
perature T tot ). We see that the solute is effectively equi- 
librated with the thermal bath of water molecules after 
2 ps. 

For such a small molecule, it is reasonable to predict 
that the impact of thermal fluctuations on the structure 
of Si5HioO(II) may have a marked impact on its optical 
properties. From previous work, [Til IT^L IT3 | we know 
that the impact of strain on Si clusters leads to large red 
shifts in the absorption gap. This is also observed in our 
finite temperature simulation. Taking snapshot struc- 
tural configurations approximately every 0.1 ps from this 
simulation, we approximate the instantaneous absorption 
gap as the difference in total energy between the ground 
state and the triplet first excited state. This is shown in 
the black curve of Fig.[^b). It is clear that variations in 
the absorption gap are highly correlated with the average 
kinetic energy (or instantaneous "temperature" ) of the Si 
atoms, shown in Fig.[S{a). The larger fluctuations in the 
absorption gap occur at higher temperatures, in the 2 ps 
equilibration phase of the simulation, where the corre- 
sponding structural fluctuations are largest, producing 
the greatest strain. Once equilibrated, the fluctuations 
in the gap reach a steady standard deviation of ~ 0.3 eV. 
The equilibrated system displays a systematic red shift 
in the absorption gap of 0.7 eV from that of the solute 
in vacuo at K. 



VIII. IMPACT OF SCREENING ON OPTICAL 
ABSORPTION OF SILICON CLUSTERS 

The observed 0.7 eV shift is not an effect due to 
polarization or dispersion from the surrounding solvent 
molecules. To verify this, we remove the water molecules 
from each structural configuration used previously and 
compute the absorption gap of the isolated Si5H 10 O(II) 
molecule. We find that the resulting unsolvated gaps [red 
curve in Fig. OJa)] follow the corresponding solvated gaps 
quite closely. We observe no systematic shift to the red or 
blue due to the presence of the solvent in the equilibrated 
phase of the simulation. This is a clear indication that 
the dominant factor influencing the absorption properties 
of the solute is the finite temperature of the system. 

It is also interesting to note that even in the highly 
non-equilibrium phase of the simulation of the first 2 ps, 



structural variations remain the dominant factor in vari- 
ations of the absorption gap. In this regime, the solvated 
and unsolvated gaps are practically the same, while fluc- 
tuating over a range of 1.5 eV. 

The absence of a systematic solvent shift, due to 
screening, for the Si cluster in water is at variance with 
the solvent shifts of polar organic solutes in polar sol- 
vents, where an experimentally observed blue shift in 
the absorption gap is expected from screening arguments. 
Recent calculations by Bernasconi et a/.0] have shown 
that a solvent blue shift for acetone in water is resolvable 
from a relatively short FPMD simulation. This contrast 
with our results is interesting and may point to funda- 
mental differences between organic and inorganic solutes 
in water. 



IX. ELECTRONIC STRUCTURE 

We analyze the electronic structure of the solvated sys- 
tem described in Sec. lVlll in order to assess the validity of 
our approach to the calculation of absorption gaps, and 
to deduce why dielectric screening does not have a large 
impact on the optical properties. Figure |B| illustrates 
the charge density associated with several Kohn-Sham 
eigenstates near the Fermi level, for an isolated Si clus- 
ter (unsolvated), for the Si cluster surrounded by water 
(solvated) and for a configuration taken from a FPMD 
liquid water simulation at 300 K (water). 

We examine first the states illustrated in Figure HJc)- 
(h). We notice that the HOMO of the solvated system 
[Figure EJe)] is localized on the Si cluster solute. This 
localization of the HOMO is typical for each of the con- 
figurations examined in Fig. 03 This is in contrast with, 
e.g., the acetone- water system examined by Bernasconi 
et al. In that work, the HOMO varied in its localiza- 
tion during the course of a FPMD simulation, at times 
reminiscent of the HOMO of the isolated acetone sys- 
tem, and at others more like the HOMO of water in the 
condensed phase. 

The LUMO of the solvated system [Figure E^f)] is en- 
tirely delocalized over the simulation cell. However, when 
we alter the occupancy of the Kohn-Sham eigenstates (as 
described in Sec. llVfl to produce the singlet or triplet first 
excited states, we notice that the HOMO* - the analogue 
of the LUMO in the occupied excited state - is localized 
on the solute. In comparison with the localized states of 
the unsolvated Si cluster shown in Figure Ela)-(d), only 
the LUMO of the solvated system differs significantly. 

Our initial motivation for considering occupied ex- 
cited state configurations was to capture the screening 
response of the solvent to variations in the solute charge 
density from its electronic ground state. We find that the 
physically different charge distributions associated with 
the LUMO and HOMO* of this solvated system are not 
significant for the various estimates of the absorption gap. 
We see from Table ITTT1 that the HOMO-LUMO estimate 
is essentially the same as the ASCF singlet estimate. 
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FIG. 6: The electronic structure of various Kohn-Sham eigenstates near the Fermi level for a specific configuration chosen at 
the 2.5 ps point from the FPMD simulation of SisHioO(II) in water outlined in Fig.0 The unsolvated systems (a)-(d) contain 
equivalent solute structures to the solvated systems (e)-(h), but have all water molecules removed. Also displayed (i)-(l) are 
states near the Fermi level for a configuration taken from an FPMD simulation of water at 300 K. All charge density isosurfaces 
(green) contain 30% of the total charge of the given state. Atom colors are as indicated in Fig. 



Our use of the triplet excited state, for reasons of ef- 
ficiency (Sec. II VI) . is justified by the close physical re- 
semblance of the charge densities of the HOMO* in each 
case. The splitting between the singlet and triplet ASCF 
absorption gaps can be associated with exchange effects 
local to the Si cluster, since the relevant gap states are 
localized on the cluster and the presence of water or the 
perturbation of the molecular structure has little effect 
on this splitting, which remains ~ 0.2 eV in all cases. 

However, the difference between the LUMO and 
HOMO* has important consequences for other physical 
quantities, such as the oscillator strength (/) for this 
transition. We find that for the HOMO-LUMO tran- 
sition, / = 0.05, while for the transition from LUMO* 
to HOMO*, / = 0.60 (refer to Fig. [^b) for definitions 
of these states). Using the LUMO* is justified given its 



almost complete (99.8%) overlap with the HOMO of the 
ground state. This order of magnitude difference in / 
will have important consequences for calculated optical 
absorption. 

The delocalization of the LUMO of the solvated sys- 
tem [Fig. E^f)] is reminiscent of the delocalized LUMO 
observed in previous DFT calculations of the electronic 
structure of liquid water. |35l |36( An example of this 
LUMO of water is provided in Fig. Efj). To interpret 
the origin of various states near the Fermi level of the 
solvated system, we compare its density of states (DOS), 
in Fig. [7| with that of the unsolvated system and that 
of pure liquid water. The energy scale in this figure is 
shifted in such a way that the HOMO of the solvated sys- 
tem is at zero. Given that the HOMO and LUMO of the 
solvated system are qualitatively similar to the HOMO 
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TABLE III: Various DFT estimates of the absorption gap 
for the Si5HioO(II) molecule: in its vacuum, zero temper- 
ature structure; at 300 K in vacuo; and at 300 K in water. 
The finite temperature estimates of the gap use the molecular 
structure of the snapshot considered in Fig. El Also shown is 
the splitting between the singlet and triplet ASCF estimates 
of the gap. All energies given in eV. 



Phase 


vacuum 


vacuum 


solvated 


Temp (K) 





300 


300 


HOMO-LUMO 


4.5 


3.7 


3.9 


ASCF singlet 


4.4 


3.8 


3.9 


ASCF triplet 


4.2 


3.6 


3.8 


singlet-triplet 


0.19 


0.18 


0.17 



of the isolated solute and the LUMO of pure water, re- 
spectively, we choose to align these pairs of states for the 
sake of comparison. It is worth noting that this seem- 
ingly arbitrary alignment is consistent for other effective 
single-particle states. For example, the lowest valence s 
states of the oxygen atoms in each system are correctly 
aligned at around ~ — 20 eV on this energy scale (this 
part of the DOS has been omitted for clarity in Fig. 0. 
We also note that the isolated LUMO of liquid water, 
shown at ~ 2.5 eV in Fig. Etc), is the subject of some 
controversy and may not reflect a physically observable 
transition. |18ll37| 

The DOS of the unsolvated Si cluster indicates another 
state quite close in energy to the LUMO of the solvated 
system, and, in fact, the LUMO+1 Kohn-Sham eigen- 
state of the solvated system resembles this unsolvated, 
localized, empty state. This is also true of the HOMO* 
of the singlet and triplet configurations. Therefore we 
conclude that occupying the first excited state leads to a 
reordering of states near the Fermi level, such that the de- 
localized LUMO is replaced by a localized HOMO*. The 
preference for a localized HOMO* is due to the binding 
energy associated with the underlying charge distribu- 
tion, which contains a "hole" (unoccupied LUMO*) lo- 
calized on the Si cluster. To further illustrate this point, 
we see that for the states near the Fermi level in liquid 
water [Fig. HDT)-(l)] there is no significant difference be- 
tween the charge densities of the LUMO and HOMO*. 
There is no localized state nearby in energy in the con- 
duction band with which to swap, and also there is no 
localized "hole" to favor such a reordering of energy lev- 
els. 



X. CONSEQUENCES FOR EXPERIMENT AND 
THEORY 

There are several conclusions that experimentalists 
may draw from our computations. With regard to reac- 
tivity, we see that for absorption processes, the existence 
of the silanone group is very unlikely, given its low ac- 
tivation energy for reaction with water. Note that the 
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FIG. 7: The density of states (DOS) of (a) the unsolvated 
Si cluster of Fig. ISfa - ); (b) the solvated Si cluster of Fig. Etc-); 
and (c) that of liquid water. The DOS of water was calcu- 
lated by thermal averaging over 20 ps of a FPMD simulation 
of 32 water molecules at ambient pressure and 300 K. The 
average broadening of states for this simulation was 0.08 eV 
which was then used to broaden the discrete spectra of (a) and 

(b) . The energy scale in (b) is shifted to place the HOMO 
at zero (vertical dashed line). The shift in energy applied to 
(a) places the HOMO at zero also, while the shift applied to 

(c) aligns the LUMO with the LUMO of (b) (vetical dashed- 
dotted line). The oxygen s-states have been omitted, but lie 
at ~ -20 eV. 



existence of a photostable silanone in the presence of wa- 
ter, proposed by Zhou and Head, |3^| is only relevant for 
emission processes. Also, other work indicates that 
the significant red shift induced by this functional group 
may be reproduced by complete passivation with other 
types of oxygen containing groups, such as hydroxyls and 
bridging oxygens. 

We expect that the large red shift caused by thermal 
strain in our finite temperature simulations (Fig. [SJ may 
be more pronounced for this small prototypical cluster 
than it would be for a larger Si QD. For larger Si QDs, 
with more rigid surface reconstructions, we expect that 
finite temperature effects will play a smaller role in their 
impact on the absorption gap. This is further reinforced 
by the small temperature dependence of the lowest direct 
gap of bulk Si, which shifts to the red by only ~ 0.05 eV 



10 



upon increasing the temperature from K to 300 K. j3^ 
The lack of any screening impact on the absorption 
of such a small Si cluster, with a relatively large dipole 
moment indicates that for larger Si QDs, this effect may 
also be negligible. The likelihood of complete passiva- 
tion of the surface of a large Si QD with oxygen will 
significantly reduce the total dipole moment of such par- 
ticles. Furthermore, the preference for oxygen contain- 
ing passivants which effectively preserve the local tetra- 
hedral symmetry around each Si atom guarantees that 
the electronic states near the Fermi level are delocalized 
over the entire QD. Therefore, screening impacts are re- 
duced, given that they are more likely to be limited to 
the surface region of the QD, not extending far into the 
core. However, in larger Si QDs a smaller red shift due to 
thermal strain may permit the observation of small sol- 
vent shifts due to screening. This issue requires further 
investigation. 

For theorists, it seems that simulating solvated Si QDs 
may not require large FPMD simulations which include 
water. Given that the screening impact is so small, there 
may be no advantage to including many solvation shells 
of water molecules in a simulation. On the other hand, 
analysis of ground state structures and active vibrational 
modes at a given finite temperature may provide more 
useful information than large finite temperature simula- 
tions and thermodynamic averaging. 

XI. CONCLUSIONS 

This is the first theoretical analysis of the impact of 
aqueous solvation on the optical properties of Si nanopar- 



ticles. While there have been many studies on the impact 
of polar solvents on the absorption properties of organic 
molecules, so far no investigation has been reported for 
inorganic solutes such as QDs. In our analysis we con- 
sidered three factors: chemical reactivity; thermal equi- 
libration; and dielectric screening. Regarding chemical 
reactivity, we find that the silanone functional group is 
extremely reactive in the presence of water and is unlikely 
to exist in aqueous solvation. We find that the bridging 
oxygen and hydroxide surface passivants are much more 
stable and therefore more probable sources of red shifts 
in the experimentally observed absorption spectra. At 
300 K, we find that thermal fluctuations induce strain in 
the small silicon cluster examined here, and we estimate 
that this strain leads to a systematic red shift of ~ 0.7 eV. 
Upon removing the screening impact of the surrounding 
water, we determine no noticeable impact of solvent po- 
larization or dispersion on the optical absorption gap of 
this silicon cluster. Given the thermal stability of larger 
QDs, we conclude that chemically stable Si QDs will have 
robust optical properties in the presence of water, further 
justifying their use as stable and efficient optical tags in 
biological sensing applications. 
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